Is a single photon's wave front observable? 
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Abstract 

The ultimate goal and the theoretical limit of weak signal detection is the ability to detect 
a single photon against a noisy background. In this situation the inescapable noise produced 
by the measuring device itself may be the main threat, but the uncertainty principle strongly 
restricts possible experimental techniques of increasing the signal-to-noise ratio. For example, a 
weak classical signal from a remote source can be distinguished from the local noise at the same 
frequency through its spatial correlations (using phase sensitive detectors; coincidence counters; 
etc) - i.e., by sensing its wave front. This method seems impossible in case of a single incoming 
photon, since it can only be absorbed one single time. Nevertheless such a conclusion would be too 
hasty. In this paper we show, that a combination of a quantum metamaterial (QMM)-based sensor 
matrix and quantum non-demolition (QND) readout of its quantum state allows, in principle, to 
detect a single photon in several points, i.e., to observe its wave front. 

Actually, there are a few possible ways of doing this, with at least one within the reach of current 
experimental techniques for the microwave range. The ability to resolve the quantum-limited signal 
from a remote source against a much stronger local noise would bring significant advantages to 
such diverse fields of activity as, e.g., microwave astronomy and missile defence. 

The key components of the proposed method are 1) the entangling interaction of the incoming 
photon with the QMM sensor matrix, which produces the spatially correlated quantum state of 
the latter, and 2) the QND readout of the collective observable (e.g., total magnetic moment), 
which characterizes this quantum state. The effects of local noise (e.g., fluctuations affecting the 
elements of the matrix) will be suppressed relative to the signal from the spatially coherent field 
of (even) a single photon (see Fig. [[]). 
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FIG. 1: The wave front of a single photon originating from a remote source being detected by 
spatially separated qubits (a quantum metamaterial-based sensor matrix) against uncorrelated 
background noise (nearby radiation sources and local fluctuations). 

INTRODUCTION 

We will illustrate the possible implementations of this scheme using a simplified model, 
an example of which is shown in Fig. [2j Here the QMM matrix is modelled by a set of N 
qubits, which are coupled to two LC circuits: the one (A) represents the input mode, and 
the other (B) the readout. This model is closest to the case of microwave signal detection 
using superconducting qubits, which is both most feasible and most interesting (at least 
from the point of view of radioastronomy) . Nevertheless our approach and conclusions 
apply generally, mutatis mutandis (e.g., to the case of doped photonic cavity QMM matrix 
in the optical range) . We will begin by discussing how such a detector system could work in 
principle. Before going on to demonstrate that a clear distinction between a single incident 
photon and the vacuum can be seen in the response of a simple two-qubit detector array 
using a fully quantum mechanical model. Then finally we will explore the role of inter- 
qubit coupling and increasing the size of the QMM array using a semi-classical mean field 
approach. 

The system of Fig. [2] can be described by the Hamiltonian 

H = H a + V a + H qb + Vb + H b + -f/noise- (1) 
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Here 



H a = u a (a^a + 1/2) + f(t)(a! + a) 



describes the input circuit, excited by the incoming field; 



H, 



qb 



1\ N 
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(2) 



(3) 



is the Hamiltonian of the qubits; 

H b = u b (tfb+ 1/2) + h(t)(tf + b) 



(4) 



is the Hamiltonian of the output circuit with the probing field, used in case of so called IMT 
readout (see, e.g., [l|); the terms 



(5) 



describe the coupling between the QMM matrix and the input and output circuits; finally, 

H noise = J2(m^ + Vj(t)a;) (6) 

3 

takes care of the ambient noise sources, which we in agreement with our assumptions take 
to be independent: (£j(t)£k(t')) oc 5jk] (£j(i)5r}k(t')) = 0). The standard way of introducing 
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FIG. 2: Schematic for the photon detector system. Photons are incident on to the QMM matrix, 
which is comprised of N qubits in this case. The QMM matrix is also coupled to the readout tank 
circuit in order to perform quantum non-demolition measurement. 
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the intrinsic noise in both LC-circuits and qubits through appropriate Lindblad operators 
in the master equation (see, e.g., can be used as easily and will be utilised in some of the 
following numerical results. The current approach allows a more transparent qualitiative 
treatment and a more straightforward accommodation for 1 / /-noise. 



A QUANTUM METAMATERIAL BASED DETECTOR ARRAY 

For a simple illustration of how a single photon can be simultaneously detected at several 
points in space, consider the case when there is one photon in the input field, which is 
coupled uniformly to two identical noiseless qubits initially in their ground states, and the 
readout circuit is switched off. In this case the system undergoes vacuum Rabi oscillations, 
and its wave function is j^j 

|tf(t)) = cos(V2</°*)|l) <8> Ui) ® I U) + 

,sin(^t)|0) 9 Iil) ^ lt2) + ' tl) ^ |i2) . (7) 

V2 

At the moments when the first term vanishes, t n = (tt/2 + ™)/V2g a , the qubits are in 
the maximally entangled Bell state, and the QND measurement of their summary "spin" 
in z-direction realizes the observation of a single photon's presence (a Fock state |1) of the 
circuit A) at two spatially separated points (locations of the qubits 1 and 2). 

A literal realization of such a scheme for observing a single photon's wavefront in multiple 
points is theoretically possible, but hardly advisable: The resonant transfer of the incoming 
photon into the qubit matrix and back is vulnerable to absorption in one of the qubits. A 
better opportunity is presented by the dispersive regime, when the mismatch between the 
qubits' and incoming photon's resonant frequencies, 5fL = \u a — */ + e 2 A ^> allows to 



use the Schrieffer- Wolff transformation to reduce the interaction term V n to 



ii.' 



Now the effect of the input field on the detector qubits is the additional phase gain propor- 
tional to the number of incoming photons, which can be read out using a QND technique. 

Let us excite the input circuit with a resonant field, f(t) = f e {t) exp[— iu a t] + c.c, with 
slow real envelope function f e (t). Neglecting for the moment the rest of the system, due to 
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the weakness of the effective coupling g 2 /SQ in (151) , we can write for the wave function of 
the input circuit 



and 



Here D(a) with 



ij t \Mt))"fe(t)(a + a^ a (t)), 



\Mt)} « e-^o W0](aW)|^ (0)) = D( a )|^ a (o)>. 



a(t) 



is the displacement operator 



D(a) 



dt'f e {t' 



^aa< —a a 



(9) 
(10) 

(11) 
(12) 



Acting on a vacuum state, it produces a coherent state, D(a)\0) = \a). Therefore, assuming 
that the input circuit was initially in the vacuum state, the average 



a*a) t ~ {ip a (t)\(T a\ip a (t)) ~ (a(t)\a* a\a(t)} = \a(t)\ 



dt'Ut' 



L-/o 



(13) 



Therefore the action of the incoming field on the qubits in the dispersive regime can be 
approximated by replacing the terms H a and V a in the Hamiltonian ([TJ with 



In the Heisenberg representation the "spin" of the jth qubit, 



(14) 



(15) 



satisfies the Bloch equations, which in case of zero bias and only z- noise (e 3 - = 0; £j(t) = 0), 
and neglecting for the moment the interaction with the readout circuit, take the form 



d 



-s^(t) = 2[ lj \a(t)\ 2 + 



dt 



S y(t) = -2[ 7 >(t)| 2 + ru(t)}s^t) - A jS *(t); 



(16) 



d 



or, introducing = ± isj 



d 



- S f(t) = T{HiM t )\ 2 + vM4( t ) + iA i sZ M ; 



(17) 
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Let us initialize the qubit in an eigenstate of <rj (i.e., in an eigenstate of unperturbed 
qubit Hamiltonian, since €j = 0). Then s z (0) = 0,s^(0) = sj(0) (i.e., 1 or -1), and the 
equations (fT7|) can be solved perturbatively: 

'(t) = exp 



Jo 



-2A^(0) 



[ lj \a{t'')\ 2 + r l] {t ,l )}dt''Yt' - 

t' 

[ lj \a(t")\ 2 + r ]j (t")}dt'dt" 



(18) 



o jo 



Assuming identical qubits identically coupled to the input circuit, we finally obtain for the 
collective variable (z-component of the total qubit "spin" of the QMM matrix) 



A? 



S*(*) = 5>*(i) « -2 7 A S *(0)iV 



/ / \ a (t")\ 2 dt'dt" + / / ±-y2 Vj (t")dt f dt" 

Jo Jo Jo Jo N - =1 



(19) 



The second term in the brackets is the result of local fluctuations affecting separate qubits 
and is therefore, in the standard way, ~ \/N times suppressed compared to the first term 
(due to the regular evolution produced by the spatially coherent input photon field). The 
variable S z can be read out by the output LC circuit, e.g., by monitoring the equilibrium cur- 



rent/voltage noise in it ^|. The signal will be proportional to the spectral density ((S Z ) Z } UJ , 
i.e. to the Fourier transform of the correlation function (S z (t + r)S z (t)). Due to the quan- 
tum regression theorem [jj, the relevant correlators satisfy the same equations (1171) as the 
operator components themselves, and the "regular" and "noisy" terms originating from ( fl9l) 
will indeed be 0(N 2 ) and O(N) respectively. 



DETECTING A SINGLE INCIDENT PHOTON 



To investigate the level of sensitivity of the proposed detector system, we consider the 
example of a QMM matrix comprised of two qubits coupled to the input mode and readout 
oscillator, as shown in Fig. [2J We assume that the input field has a given number of photons 
incident on it and is initially found in a coherent state, |a), with an average of | or | 2 photons 
and therefore take f(t) = 0. We also take h(t) = and assume that the intrinsic noise in 
the detector system is sufficient to drive the readout field and allow detection of the incident 
photons. 
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In order to fully account for the effects of decoherence and measurement, we make use of 



the quantum state diffusion formalism [6| to describe the evolution of the state vector 



\di/>) = -iH \vj>) dt + 



dt+ 

%> i^)de, (20) 



where \dijj) and dt are the state vector and time increments respectively, Lj are the Lindblad 
operators and d£ are the stochastic Wiener increments which satisfy d£ 2 = d£ = and 



d£d£* = dt. In this case we take 

H = H a + V a + H qb + V b + H b . (21) 

To replace H no i se and model the natural effects of decoherence on the qubits we instead 
have the Lindblad operators L z = y/2T z (T^} and L xy = ^/2T xy a^cr_ acting on both qubits. 
These operators describe relaxation in the ^-direction and dephasing in the x-y plane of 
the Bloch sphere respectively. To account for the weak continuous measurement of the 
output field we also take L b = ^j2T b b. From the real and imaginary parts of (b) we can 



extract the expectation values for the position x b = y/T/2uJ~ b (b + tf) and momentum p b = 
i^Ju b /2 — b) operators. 

We solve flU numerically with e = u a = 1, A = 0.1, g] = 0.1, g) = 0.01, T z = T xy = 10~ 4 
and T b = 0.1. To ensure that we are in the dispersive regime we take u b = 2.5 to give 
SQj = 1.495 3> g a y For the input and output field Hamiltonians we truncate the Hilbert 
space to the lowest 30 states. We assume that the output field begins in the vacuum state 
and both qubits are initially in an eigenstate of aj, (|0) + |1)) /y/2. We investigate the cases 
where the input field begins in the vacuum state and coherent states with an average of 1 
and 2 photons. In each case the expectation value is integrated over 1000 periods of u a and 
then the Fourier transform is taken to produce the power spectral density for the position 
and momentum quadratures. 

An example of the power spectral densities for the x and p quadratures of a typical 
experimental trajectory are shown in Fig. [3j We can see that in both cases there is a clear 
distinction in the readout depending on whether a photon is incident upon the detector or 
not. When the input field begins in a coherent state the peak in the power spectral density 
is shifted to lower frequencies and decreases in magnitude compared to the clear sharp peak 
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FIG. 3: Power spectral density for the position (^(xi ) )' 2 J and momentum y(j>fe) y quadratures of 
the readout tank circuit for the cases where the detector has 0, 1 and 2 photons incident upon it. 

seen when the input field is in the vacuum state. The response is smeared out across the low 
frequency region leading to a higher average power when there photons there are photons 
incident upon the detector. This is particularly clear in the case of the position operator 
where there is a clear distinction of approximately one order of magnitude. 
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SCALING OF THE QUANTUM METAMATERIAL SENSOR ARRAY 



In order to further investigate the role of an increasing number of qubits and of interqubit 
couplings in the QMM matrix, we now consider the following reduced Hamiltonian 

H = -\ E + + a E °i ( 22 ) 

with the qubits driven by a common harmonic off-resonance signal (modeling the input V a 
of Eq.(PP)) and local noise coupled through a z : 

6j {t) = e sm(ut) + V2D^ (t) . (23) 

Here (£(t)) = and (£j(t)£i(t') = 5ji5(t—t'). This treatment is consistent with our qualitative 
approach of Eqs.flMHJ)- 

In the case of N uncoupled qubits, we can describe the system by N independent master 
equations, 

(24) 

for a single-qubit density matrix, and average the observable quantities. These results are 
shown in Fig. H^,,b. The spectral density amplitude of of the ^-component of the total 
"spin" (i.e., square root of the spectral density ((S z ) 2 ) u ) demonstrates a small, but distinct 
peak due to the external drive, in addition to the large noise-driven signal. The increase 
of the number of qubits, predictably, increases the signal to noise ratio. The increase is 
in qualitative agreement with the \/N behaviour, though numerically somewhat smaller 
(approximately doubling rather than tripling as iV increases from 1 to 9; see Fig. [4^, inset). 

The introduction of qubit-qubit coupling also increases the signal to noise ratio. In this 
case, we solve the master equation ( )24|) for two coupled qubits, using the generalized Bloch 
parametrization of the two-qubit density matrix: 

P=\ E n «* a « ® a b- (25) 

a,b=0,x,y,z 

The results in Fig. |4fc show that while the overall signal amplitude is suppressed by qubit- 
qubit coupling, the relative amplitude of the signal significantly increases. 



ap 
dt, 



H(t)J 
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FIG. 4: (a) Spectral density of total detector "spin" S z (the square root of the spectral density 
of fluctuations, ((S Z ) 2 ) UJ ) in a qubit in the presence of noise and drive. The signal due to drive 
is a small thin peak on the left of the resonant noise response. Inset: Signal to noise ratio as the 
function of number of qubits. (b) Same in case of 8 qubits. Inset: A close up of the signal-induced 
feature. The noise is suppressed in case of 8 qubits (blue) compared to the case of a single qubit 
(red), (c) The spectral density of S z in case of two coupled qubits. Note that the significant shift of 
the resonant frequency of the system (position of the noise- induced feature). Inset: Signal response 
amplitude (left) and signal to noise ratio (right) as functions of the coupling strength. 



11 



DISCUSSION 



Though the possibility to observe a single photon's wave-front requires the detection of 
a weak, remote signal against the background of local fluctuations, the standard signal- 
to- noise ratio v^V-enhancement due to the A-element coherent uncoupled QMM matrix is 
unlikely to be of much practical use. Noticing that the effect of the input field is nothing but a 
simple one-qubit quantum gate applied to each element of the matrix and that introducing a 
simple qubit-qubit coupling scheme can improve matters, we can ponder a more sophisticated 
approach. By performing on a group of qubits a set of quantum manipulations, which would 
realize a quantum error correction routine, one can hope to improve the sensitivity of the 
system. We will consider this approach in a separate paper. Another possibility is to consider 
different types of unit elements in the QMM array. For instance, an array of SQUID rings 
would offer the potential for significant frequency conversion between the incoming source 
and measurement circuit by either up or down conversion Q. 

In conclusion, we have shown the possibility in principle to detect the wavefront of a 
single photon using the quantum coherent set of spatially separated qubits (a quantum 
metamaterial sensor matrix). The key feature of this approach is the combination of the 
nonlocal photon interaction with the collective observable of the QMM matrix and its QND 
measurement. Besides the intriguing possibility to test the limits of application of quantum 
mechanics, the realization of our approach would allow to greatly improve the sensitivity of 
radiation detectors by suppressing the effects of local noise as well as lowering the detection 
barrier to the minimum allowed by the uncertainty principle. 
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